


m 



Imaging the Earth's Interior: the Angular Distribution of Terrestrial Neutrinos 

Brian D. Fields-f^'^ and Kathrin A. Hochmuth*^ 

^ Center for Theoretical Astrophysics, Department of Astronomy, University of Illinois 
^Department of Physics, University of Illinois 



Abstract 
Decays of radionuclides throughout the Earth's interior produce geothermal heat, but also are a source of 
antineutrinos; these geoneutrinos are now becoming observable in experiments such as KamLAND. Hereto- 



o 
o, 

Cn I fore, the (angle-integrated) geoneutrino flux has been shown to provide a unique probe of geothermal heating 



due to decays, and an integral constraint on the distribution of radionuclides in the Earth. In this paper, 
we calculate the angular distribution of geoneutrinos, which opens a window on the differential radial distri- 
bution of terrestrial radionuclides. We develop the general formalism for the neutrino angular distribution. 
We also present the inverse transformation which recovers the terrestrial radioisotope distribution given a 
measurement of the neutrino angular distribution. Thus, geoneutrinos not only allow a means to image the 
_ Earth's interior, but offering a direct measure of the radioactive Earth, both (1) revealing the Earth's inner 

^^ ' structure as probed by radionuclides, and (2) allowing for a complete determination of the radioactive heat 

f^ , generation as a function of radius. Turning to specific models, we emphasize the very useful approximation 

^— N , in which the Earth is modeled as a series of shells of uniform density. Using this multishell approximation, 

r^ ' we present the geoneutrino angular distribution for the favored Earth model which has been used to calculate 

I ' geoneutrino flux. In this model the neutrino generation is dominated by decays of potassium, uranium, and 

Ph! . . 

^ thorium in the Earth's mantle and crust; this leads to a very "peripheral" angular distribution, in which 2/3 

of the neutrinos come from angles 9 > 60° away from the downward vertical. We note that a measurement 
. . - of the neutrino intensity in peripheral directions leads to a strong lower limit to the central intensity. We 

%^ ' also note that there is some controversy about the abundance of potassium in the Earth's core; different 

geophysical predictions lead to strongly varying, and hence distinguishable, central intensities {0 < 30° from 
the downward vertical). Other uncertainties in the models, and prospects for observation of the geoneutrino 
angular distribution, are briefly discussed. We conclude by urging the development and construction of 
antineutrino experiments with angular sensitivity. 
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I. INTRODUCTION 

The decays of radioactive species within the Earth generate an important component of geother- 
mal heat. However, a quantitative accounting of the radioactive energy generation of the Earth 
requires a detailed knowledge of the abundance distribution of the key long-lived radioisotopes — 
uranium, thorium and potassium-inside the Earth. Because our knowledge of these abundance 
distributions is incomplete, the radiogenic heat output is consequently model-dependent and thus 
uncertain. A fundamental diagnostic is the "Urey ratio" which measures the ratio -PradZ-Piost of 
total radioactive heat production to the surface heat loss. Current estimates span the range [1], [2] 

Pr.d/P\ost ~ 0.5 - 0.6 (1) 

which seems to suggest that radiogenic heating is a dominant heat source, but not the only one. It is 
even possible that the Urey ratio is closer to 1, as the primordial heat from Earth's formation should 
have radiated away a long time ago [3], and so the Earth's heat should be fully radiogenic. This 
would imply that more radioactive material is hidden in the Earth. The strength of these qualitative 
conclusions thus hangs on the strength of the quantitative measurements of the radiogenic heat 
production and total heat loss. 

A beautiful if challenging means of measuring the radiogenic heat production of the Earth 
follows from the realization that /3-decays not only are a heat source but also produce (anti)- 
neutrinos. Decades ago, the prescient work of Eder [4] and later Krauss, Glashow, and Schramm 
[5] pointed out that the radioactive heat flux and the neutrino flux from the Earth are tightly 
linked. A measurement of the neutrino flux would constrain the radiogenic heat production, and 
thus offer a new and direct measure of the Urey ratio. Clearly, a precise measurement of the Urey 
ratio would provide important insight into the interior structure and dynamics of the Earth. 

This vision of neutrino geophysics has enjoyed a major advance recently, with the first detection 
of geoneutrinos by the KamLAND Collaboration [6]. In fact, the detection of these geoneutrinos 
was only a side product of the KamLAND detector, as it is primarily designed to detect the flavor 
change of reactor antineutrinos produced by the Japanese nuclear power plants. The KamLAND 
Collaboration reported a total of 9 geoneutrino events, of which they estimate 4 are from 238 ^ 
decay, and 5 are from ^^^Th decay [6]. Obviously these are early days, but we are encouraged that, 
even in the prescence of a dominant anthropogenic background, KamLAND has demonstrated that 
geoneutrinos exist at observable levels. 

Several groups have taken interest in these neutrinos and have in conjunction with geophysi- 



cists and geologists worked on models for the distribution of radioactive elements in the Earth 
and predicted the (angle-integrated) geoneutrino flux at different detector sites [7-11]. The flux 
calculations have already become rather sophisticated, and are based on detailed geophysical mod- 
els, in some cases even including the effects of the anisotropic radioisotope densities in the crust. 
These models confirm in detail that indeed the geoneutrino heat flux is proportional to radiogenic 
heat flux, but with the important caveat that the exact scaling between the two depends on the 
detailed abundance and density distributions within the Earth. ^ In this sense, the geoneutrino flux 
measurement is an integral measure of the radioisotope distribution. 

In this paper we want to consider the angular distribution of the geoneutrino flux. We show that, 
once well-measured, the angular distribution can be inverted to recover the full density distribution 
of radionuclides-a tomography of the structure and radiogenic heat generation of the Earth. In 
addition, we come to the conclusion that with a future low-energy antineutrino detector with even 
crude angular resolution, it will be possible to distinguish between the different Earth Models and 
solve the problem as to how much radioactive material is contained in the Earth and where it is 
located. Thus, the angular distribution provides a differential measurement of the radioisotope 
distribution, and can reveal a wealth of new information about the structure and content of the 
Earth. 

We first present the formal calculations in §11. After outlining the general formalism, we consider 
the useful approximation of a uniform density shell, from which a multiple-shell model of the 
Earth can be constructed. We then review (§111) models for the radioisotopic content of the Earth. 
Using these models we construct physically motivated illustrations of plausible geoneutrino angular 
distributions (§IV). We present several plots for different geological abundance predictions, finding 
in general that the high radioistope content of the crust leads to a large "peripheral" geoneutrino 
signal. We investigate the possibility of ^''K in the core, as recently suggested by several groups 
[12-14]; if core ^'^K is at the high end of these predictions, then the resulting neutrino signal 
could be quite large and should lead to a readily observable central intensity peak. Finally, we 
present conclusions §V, and discuss the exciting possibilities that will arise when we are able to 
use neutrinos to image the interior of the Earth. 

Before we begin the formal development, a word of clarification seems in order. Note that for 



^ It is true that, at least in a spherically symmetric Earth, the surface heat flux is directly related, via Gauss' law, to 
the total mass of radionuclides (c.f. eq. 24). However, a simple Gauss' law argument fails in the case of neutrinos, 
because at each point they are emitted isotropically, not just radially, and as we will see, the "sideways" neutrinos 
make a large contribution to the total surface geoneutrino flux. 



brevity, we will refer to the emitted particles as neutrinos, although they are of course antineutrinos. 
In fact, there is some regular z^g production due to '^'^K electron captures. However, the branching 
here is 10.72 % of all ^'^K decays, and thus the majority of the ^''K nuclei /?-decay and yield Pg with 
a continuous energy spectrum. We look forward to the day when the monoenergetic ^''K electron 
capture Ve flux (and angular distribution!) can be measured and compared against the i>e signal. 
However, in this paper we will consider only the dominant, /3-decay i>e signals from K, Th, and U. 

II. FORMALISM 

The fundamental quantity we wish to calculate is the differential intensity I, or surface bright- 
ness, of geoneutrinos at the surface of the Earth. This is just the distribution of neutrino flux F 
versus solid angle: I{9,(j)) = dF/dD.. Both angular coordinates are local and observer-centered: 
6 G [0, 7r/2] is the nadir angle, i.e., the angle measured from the downward vertical (so that the 
center of the Earth is at = 0, and the horizontal is at = 7r/2; see Figure 1). The angle (f) is 
an azimuth. These angles thus cover the "sky" underfoot (or rather, the terrestrial hemisphere) 
which we wish to image. In this paper we will consider only the case of a spherically symmetric 
Earth (including the outermost layers). This guarantees that the intensity / = I {9) is azimuthally 
symmetric and so only depends on the nadir angle. In this case, we have I{6) = {2TT)^^dF/dcos6. 

o 




FIG. 1: The basic geometry of the problem. An observer at point O measures an intensity, which sums the 
emission along a line of sight s. The nadir R^ points to the center C of the Earth, which makes an angle 9 
with the line of sight. A given point along s is at geocentric distance r. 



The neutrino intensity in any direction depends on the distribution of sources along that line of 
sight. The governing equation is that of radiation transfer for neutrinos, which is formally identical 
to the usual expression for photons [15], although of course the microsphysics is quite different (the 



Earth is optically thin to neutrinos, but neutrinos do undergo oscillations). Then over a line of 
sight s, ignoring scattering and absorption, the intensity changes according to 

dl/ds = q{s)/47r (2) 

where q is the source function at point s, which measures the local neutrino production rate per 
unit volume. 

For each radioisotope species i, this takes the form 

qi = — = • (3) 

Ti TiUli 

Here rii and pi are the local number and mass densities, respectively, of species i in the Earth. The 
Ti is the mean lifetime and rrii is the mass of an i nucleus. The total source is just a superposition 
of all species: q = ^i qi- The effects of neutrino oscillations are not yet included; we will address 
this below (§11 B). 

Integrating eq. (2) over a line of sight at nadir angle 9, we find the intensity 

r2R^ cos 9 

m= q{s)ds (4) 

JO 

where sis centered on the observer, as seen in Figure 1. 

Since models of the Earth's structure and composition are expressed in terms of the radius r, 
it is very useful to rewrite eq. (4) in these geocentric coordinates. From Figure 1 we see that the 
s = r+ R® , so that f = s— R^ and thus 



s{r, 9) = iJ© cos ± •y/r2 - (1 - cos2 9)R% (5) 

where it corresponds to the far side (near side) of the midpoint s = R^ cos 9 of the line of sight. 
Therefore, for a fixed nadir angle 9, we can transform the integral to the geocentric coordinate 
system: 

/.(«)= 2/"- 4r , "-W- (6) 

J R^ Sine J r^ - {I - cos^ 9)Rl 



where qi{r) can in a non-spherically symmetric case also be a function of (geocentric) latitude and 
longitude. In this case, the factor of 2 is replaced by the sum of integrals for the near-side and the 
far-side of the midpoint s = Rcos9. In our special case of spherical symmetry, the contribution 
from the near half of the sphere is the same as the contribution from the far side. 

Equation (6) explicitly demonstrates that the intensity I{9) is an integral transformation of 
the source distribution q{r). In fact, this mapping is a form of the Abel transform [16], which is 



used in deprojection problems in both astrophysics [17] and geophysics [18]. Thus it is clear that a 
determination of the intensity distribution offers a measure of the source distribution. Namely, one 
can invert the transformation (deproject the image) to fully recover the complete source distribution 
q{r); this inversion procedure is explicitly presented in Appendix A. In other words, a measurement 
of the angular distribution of geoneutrinos not only yields an image of the Earth's radioactive 
interior, but this image can also be inverted to give a tomography of the terrestrial radioisotope 
distribution. Clearly, the geoneutrino angular distribution offers a unique and powerful probe of 
the interior of the Earth. This power will be further illustrated below, where we will see that even 
a partial (low-resolution) determination of I{9) offers important geophysical information. 

In evaluating eq. (6), it will be convenient to introduce dimensionless scaled variables: a radial 
fraction x = r/R^ G [0, 1], a local mass fraction Oj = Pi/p, and a local density measure p = p/p 
(with p = SM/AirR^ the mean Earth density). It will also be useful to denote the nadir angle 
cosine as /i = cos^. We then rewrite eq. (6) as a product of two terms 

I^{e) = Ii,og^{p) (7) 

Here the dimensionful overall magnitude is set by 

where Ni is neutrino multiplicity, i.e., the number of geoneutrinos released per decay chain. Values 
of li^Q appear in Table I. The dimensionless angular distribution or (akin to the "phase function" 
of radiation transfer [15]) is the heart of this paper, and is contained in the function 

m{p)= f'dx f(;]" (9) 

J^/^^ \JX^ - (1 - p^) 
which requires knowledge of the density distribution of each radioisotope, usually presented in 
terms of mass fractions Oj and a total density profile p via pi = Pi/pi = {ai/di)p. 

The flux of geoneutrinos integrated over different annuli is also of interest. We quantify this in 
terms of the flux exterior to the nadir angle 9, via 

Fi{> 9) = Fi{< p) = f d^h{p,4,) (10) 

= 27r/i,o r dpg{p) (11) 

Jo 

= 27Tl,^oH{p) (12) 

where second and third expressions assume spherical symmetry, and where the dimensionless quan- 
tity H encodes the angular dependence. A consequence of this definition is that the total neutrino 



flux is given by Fi{9 > 0) = Fi{fi < 1) 
calculations [5, 8, 10]. 



2TrIi^QH(l). This can be compared with existing 





TABLE I: Properties of the Princ 


iple Geoneutrino 


Source Nuclei 




Radioisotope 


mean life 


i/e multiplicity 


isotopic 


mean terrestrial 


intensity normalization 




species 


r (Gyr) 


N, 


abundance 


abundance di 


lifi [neutrinos cm^^ s"^ sr" 


-'] 


40k 


1.84 


1 


0.0117% 


1.8 X 10"® 


2.4 X 10^ 




238U 


6.45 


6 


99.2745% 


5.3 X 10-® 


0.48 X 10^ 




232Th 


20.3 


4 


100% 


1.35 X 10-® 


0.56 X 10^ 





A. The Uniform Shell Approximation 



We want to consider a single shell with a density pi^ in species i which is constant between 
Hn = 3;in-R© and rout = a^out-Re; and zero otherwise: 



Piix) 



, otherwise 



(13) 



This form allows us to simplify the integral and solve it analytically. In particular, owing to the 
constant density, the intensity at each line of sight is just proportional to the shell path length 
/\s{6) along that sightline. 




< M < Min 



FIG. 2: Geometry for the uniform shell model. The inner and outer shell radii arc r-m and rout, respectively. 
The lines of sight tangent to the inner and outer edges of the shell arc shown, and are at nadir angle cosines 
/Xin and ^out , respectively. These lines of sight divide the terrestrial hemisphere into three distinct regions, 
which are labeled. 



There are three cases to be distinguished, as one can see in Figure 2: 



aifj-) 



Pifi I VA*^ - ^out - \//^^ - M?„ I , fJ-y fJ-in 



0, 



/"out < ^i < fJ-i: 
/U < /"out 



(14) 



where 



A^i: 



1 



2 



1 _ ^2 



(15) 



These are integrable, and give 

^ [^out KP-I Pont) - flf^h{fl/fli^)] , p> Pin 

^pl^t Hp/Pout) , Pont < P< Pin 

^ < Pont 



H{p) 



(16) 



where 



-| -I 

h{u) = —u\Ju^ — 1 hi ( n + V M^ — 1 



(17) 



Note that eq. (16) cahs h{u) only in the domain n > 1, and that /i(l) = 0. The total flux for a 
uniform shell was first calculated in [5], and one can easily show that F(/i < 1) reproduces their 
result. 

There are four cases we want to illustrate; these are plotted in Figure 3. In the "uniform 
Earth" model, the density is the same throughout the whole planet (xjn = 0, Xout = 1)- This 
gives an angular distribution which grows linearly with /i, g(^p) = Pi.oP, and hence I{0) oc cos(^). 
The integrated flux which increases quadratically as H{fi) = pififi'^/2, i.e., H(6) oc cos'^9. This 
therefore gives a very centrally bright distribution. In the "uniform core" model, the density is 
non-zero only in a central region which extends from x-m = to Xout- As seen in Figure 3, this gives 
an inner intensity distribution which is similar to the uniform Earth model, as one would expect, 
but which goes to zero, as it should, at the outer tangent ^out- 

As we will see below, it turns out that a more physically motivated case is the "uniform crust" 
model, where only the layers of the Earth (from Xin to Xout = 1) contribute. This yields an intensity 
distribution 



9{p) 
H{p) 



Pifi 



P \ P 



Ph 



, P> Pin 
PifiP , P< Pin 

^ [pint Hp/Pont) - Pfn Hp/Pin)] , P > Pir. 
^P^ , P<Pir. 



(18) 
(19) 
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FIG. 3: The angular distribution of intensity for different cases of a 1-zone, uniform density model, shown 
as a function of nadir angle 6. The units in the vertical scale are arbitrary, and the normalizations of the 
different curves are chosen for clarity. 

which is also linear in /i, and thus scales as / oc cos(0), for large angles. The intensity peaks at 
inner tangent point /iin, where the line of sight is longest. A measurement of this peak would thus 
give the position of the inner edge. Interior to /^in, the intensity drops to a minimum at the nadir, 
where the column density is the smallest. But the central intensity is nonzero, and is simply related 
to the peak intensity via 



/peak _ -f(/^in) _ Ag(^i,: 



'l + Xi, 



(20) 



-^center -^(1) Argholl V ^ ~ ^in 

This useful relation allows one to use the peripheral flux, due to emission from Earth's outer layers, 
to put a lower limit on the contribution of these layers to the central flux. Any observed central 
flux in excess of this limit must be due to emission from the inner Earth. 

The last model illustrated in Figure 3 is a single shell with arbitrary inner and outer radii. The 
intensity distribution is an amalgam of the features seen in the special cases of the crust and core 
models. As with the crust model, the intensity peaks at the inner tangent point, where the line of 
sight is longest. As with the core model, the intensity goes to zero at the outer tangent. Thus, a 
measurements of the peak would give the position of the shell's inner edge, while a measurement 
of the cutoff would give the outer edge. 



B. Towards a Realistic multishell Model 

Based on the idealized calculations in the last section we want to build now a model that is 
applicable for the Earth's Interior. The Earth of course does not have a constant density. However, 
we can approximate the density structure as a series of uniform density shells. Indeed, Earth models 
are in practice typically expressed in this manner. Formally, the generalization is trivial, thanks 
to the lack of neutrino absorption and scattering, which guarantees that superposition holds: 

hie) = Y. 4'V) (21) 

shells j 

F,{>e) = ^ F^^'\>e) (22) 

shells j 

where l\ {6) and F^ {6) are the intensity and flux, respectively, from shell j. 

Moreover we must consider the effect of neutrino oscillations [9, 11]. As the mass eigenstates 
of the neutrinos differ from the flavor eigenstates, some of the initial i^e signal will be converted 
into other flavors during propagation, and thus reduce the i>e signal in the detector. In general, 
the oscillations will have a "vacuum" contribution, but will be modified due to the MSW effect, 
which is density- and energy-dependent. If the MSW effect is important, then it complicates 
the determination of the angular distribution, introducing a density- and path-dependent MSW 
suppression factor into eq. (6) and its descendants. This factor would complicate the inversion of 
the angular distribution I{6) to recover the radioisotope source distribution q{r). 

However, geoneutrinos have a very low energy and as the Earth's density is much lower compared 
to stellar objects we do not expect, that the MSW effect plays any greater role. This can be justified 
if the vacuum oscillation length L^ is much smaller than the oscillation length in matter Lm- If we 
take the values of ref. [19] for the matter in the Earth and a value for A?72^ ~ 10"^ eV^, we obtain 
Lv ~ 10 km and Lm ~ 1000 km. Therefore \L^\ <C |im|) and we see that, roughly speaking, the 
vacuum oscillation length is short compared to both the the matter oscillation length and to the 
changes in density; thus the vacuum oscillations will wash out any matter effect and average out 
the pathlength dependence. We thus follow Mantovani et al. [10] and introduce only a density- 
independent oscillation probability of (1 — 2sin^2??i2) = 0.59 in our equations. Here •d\2 is the 
electron-neutrino v^ — Vx mixing angle; solar neutrino experiments and KamLAND are best fit [20] 
by sin2t?i2 = 0.315 ± 0.035 (and hence sin^ 26'i2 = 0.81). and Am^ = (7.3 ± 0.8) x 10^^ eV^. 
Therefore the intensity scaling of eq. (8) changes to: 

^^o = (l-^sin2 2^12)40 . (23) 
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Table I sums up the important properties of the radioactive elements in question. 
In the next section we want to use this and discuss different Earth Models. 

III. EARTH MODELS 

The geoneutrino intensity depends on the radioisotope density structure pi{r), which is usually 
presented as abundances aj and a total density p, where pi = aip. The Earth's interior structure 
and total density are primarily probed via the propagation of seismic waves. These results have 
been synthesized by Dziewonski and Anderson in their Preliminary Reference Earth Model [21]. 
Following [10] we will adopt this model. In addition, for our spherically symmetric study we took 
the last three kilometers of the Earth to be sediments. 

The distribution of radioisotope abundances Oj are obtained with different geological measure- 
ments which, for most of the Earth's interior, are necessarily indirect. Consequently, the abundance 
distribution remains model-dependent. Indeed, the measurement of the angular distribution of 
geoneutrinos (as well as the geoneutrino energy spectrum) would provide a powerful new method 
to measure the radioisotope distribution. For the purposes of illustration here, we will adopt the 
values of a^ given in the reference model of Mantovani et al. [10], which are very detailed and 
draw on a wealth of geophysical data. This model takes into account the bulk silicate Earth 
model, which describes the composition of the crust-mantle system. Table II shows the adopted 
abundance distribution. Note that these are elemental abundances; we assume that the isotopic 
fractions of Table I hold throughout the Earth. This correction is particularly important for ^'^K. 

TABLE II: Mantle and crust elemental abundance distribution [10]. 



Region 


radii [km] 


a(U) 


a(Tli) 


a(K) 


Lower Mantle 


3480-5600 


13.2 X 10-3 


52 X 10-9 


1.6 X 10-"^ 


Upper Mantle 


5600-6291 


6.5 X 10-9 


17.3 X 10-9 


0.78 X lO-'' 


Oceanic Crust 


6291-6368 


0.1 X 10-^ 


0.22 X 10-6 


0.125 X 10-2 


Lower Crust 


6291-6346.6 


0.62 X 10-s 


3.7 X 10-6 


0.72 X 10-2 


Middle Crust 


6346.6-6356 


1.6 X 10-s 


6.1 X 10-6 


1.67 X 10-2 


Upper Crust 


6356-6381 


2.5 X 10-6 


9.8 X 10-6 


2.57 X 10-2 


Sediments 


6368-6371 


1.68 X 10-6 


6.9 X 10-6 


1.7 X 10-2 


Oceans 


6368-6371 


3.2 X 10-9 





4.0 X 10-4 



The composition of the Earth's core deserves special attention. The Earth's core consists largely 
of iron. But its density is lower than one would expect if the core were pure iron. Therefore it is 
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assumed that light elements in the form of alloys are present [22]. So far it was generally believed, 
that the core does not hold any significant amount of radioactive elements, as there was no evidence 
that the radioactive isotopes in question could form alloys with iron. Hence, the Mantovani et al. 
[10] reference model places radioactive elements only in the mantle and in the crust, and are absent 
from the core. 

On the other hand it is puzzling why solar chondrites have a K/U ratio that is eight times 
higher than in the crust-mantle system [23]. In this connection, it is noteworthy that recent 
experiments demonstrate that potassium does form alloys with iron under high temperature and 
pressure conditions which likely were present at Earth's formation. The maximum possible amounts 
of potassium in Earth's core suggested by the experiments range from 60-130 ppm [12], to 1200 
ppm [14], to as high as 7000 ppm [13]. In light of these experiments, we will consider the possibility 
of an additional radiogenic contribution from the core, and quantify the impact of core ^''K on the 
reference model. 

The amount of radioactive material in the core contributes also to the radiogenic heat of the 
Earth. Presently the Earth's surface loses about 40 TW or 87 mW/m^. The amount of surface 
radiation depends on convection and conduction properties of the Earth's interior. The present 
day heat production H in units of TW of uranium, thorium and potassium with a total mass in 
kg of M\], Mxh, Mk respectively is 

// = V 6.M, ^ 10 (^^\ + 2.7 (^^\ + 3.4 X 10-^ {^^] TW (24) 

where e, = Qi/rriiTp^i is the specific non-neutrino energy loss per nucleus, and where Mk is the 
total potassium mass, for which the ^'^K isotopic fraction appears in Table I. To obtain the heat 
production of potassium we used, that in 89.28% of all cases a ^''K nucleus /^-decays with a average 
energy of 0.598 MeV [24]. The other decay mode of potassium is electron capture with an energy 
of 1.505 MeV. 

With knowledge of the radioisotopic content of the Earth from geoneutrinos, eq. (24) can be 
compared with the global heat flux, and used to determine the Urey ratio (eq. 1). While the 
radioisotope abundances remain uncertain, one can rather turn the problem around, and use the 
global heat flux with eq. (24) to set an upper limit to the radioisotopic content and thus to the 
geoneutrino flux. Such an analysis has been carried out by Mantovani et al., who find that this 
"maximal radiogenic" model leads to fluxes about twice the level of their reference model. 

In the next section we will discuss the plots we created for the different Earth models. 
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IV. RESULTS 

We now combine our general formalism with various Earth models to arrive at predictions for 
the geoneutrino angular distribution. We first consider the reference model, then its variants and its 
uncertainties, and finally we comment on the effect of anisotropics in the radioisotope distributions. 

A. The Reference Model 

The reference model of Mantovani et al. [10] serves as our standard and fiducial case. The 
geoneutrino intensity distribution for this model appears in Figure 4. We see that the total intensity 
is peaked near the horizon, at large nadir angles. The strikingly "peripheral" character of this 
neutrino distribution is a direct consequence of the location of the radioisotopes in the mantle and 
crust. 

The experimental ability to detect this pattern is perhaps best quantified in Figure 5, which 
displays the cumulative angle-integrated flux F for the reference model (c.f. eq. 10). The change 
in the normalized flux F/Ftot over any angle interval gives the contribution of that interval to the 
total flux. Figure 5 thus shows that fully 2/3 of the total flux arrives in the outermost nadir angles 

e > 60°. 

This model can thus be easily tested with an experiment having even the poorest angular 
resolution. For example, an experiment with 30° resolution could test whether the counts in the 
outer 6 > 60° are a factor ~ 2 higher than the counts in the inner 6 < 60° . The confirmation of this 
model would vindicate the idea that U, Th, and K are indeed lithophilic elements that congregate 
in the mantle and crust. 

Figure 4 also shows the expected angular distribution separately for uranium, thorium, potas- 
sium and the cumulative intensity for the reference model. The uranium to thorium ratio stays 
approximately the same, whereas the amount of potassium increases steeply in the crust. The 
intensity shows a double peak, which is due to the fact, that the abundance of radioactive elements 
is lower in the middle crust than in the upper and lower crust. 

We note that the reference model gives a geoneutrino distribution qualitatively similar to the 
uniform crust model presented in §11. This of course traces to the positioning of the radioisotopes 
in the outer Earth. Of course the true geoneutrino behavior is more complicated, but it does 
appear that for a crude analysis (or in the presence of crude data) the uniform crust model would 
be a useful analytic approximation. 



13 



10 






B 

u 

m 
O 






uranium 
thorium 
potassium 
total tlux 




— r~ 
20 



40 60 

Q (degrees) 



FIG. 4: The figure shows the reference model of ref. [10]. The dashed-dotted Ime represents the total 
expected intensity. No contribution coming from the core has been added yet. The other curves show the 
intensities from thorium, uranium and potassium separately. It can be seen that the major contribution is 
coming from potassium. 

Finally, we wish to comment on the effects of sediments. We find that if the amount of sediments 
is increased by a factor of 10 the overall shape of the spectrum does not change except for a sharp 
peak towards 90°. An increase by a factor of two only increases the second peak slightly. Therefore 
we conclude, that the location of the detector will not change its ability to detect the structure of 
the mantle and crust. 

B. Models with Core Potassium 

We now turn to the possibility that the Earth's core might contain significant amounts of ^^K. 
Based on the experimental evidence that potassium can form alloys with iron we use the fiducial 
values obtained by different authors and add them to the reference model. We are aware that if 
the abundance ratios and mass ratios of the reference models are correct and only the total mass 
estimates of U, Th and K are incorrect, then the intensity of the neutrinos coming from the crust- 
mantle system should go down, but the overall shape should stay the same. This decrease in the 
peripheral intensity will make the core contribution even more dominant. On the other hand, it 
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FIG. 5: Cumulative flux expected for the reference model. The flux is normalized to the total, and thus is 
dimensionless, spanning the range to 1. 

is possible that the assumptions of the reference model are correct (incorporating the bulk silicate 
Earth model with a heat production of ~ 20 TW), but the Urey ratio is closer to 1 with a large 
amount of potassium in the core. In this case it is legitimate to add the core contribution to the 
reference model. In case of a fully radiogenic model with 20 TW coming from the crust-mantle 
system one would expect a total of ~ 330 ppm of potassium in the core. Figure 6 shows the angular 
distribution for different amounts of potassium in the core as could be found in the literature, while 
Figure 7 shows the cumulative flux for the same models. 

From Figures 6 and 7 it is clear that the introduction of '^'^K in the core can significantly alter 
the geoneutrino angular distribution. The effect is to enhance the central intensity {9 < 30°), 
possibly also raising the overall detected flux. The departure from the reference model depends of 
course on the core potassium abundance. The value given in [13] of 7000 ppm potassium in the 
core is only an upper limit, but in this case the core would clearly dominate the distribution. For 
1200 ppm of potassium in the core [14] the maximum intensities coming from crust and core are 
approximately the same. But even the much lower value obtained in [12] of 60-130 ppm might still 
be detectable with a future neutrino detector. 

A future low energy neutrino detector with angular resolution will be able to distinguish be- 
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FIG. 6: The curves show the possible abundances for potassium in the core as found in [12]. [13], [14]. The 
dashed-dotted line is the intensity of the reference model. 

tween the opposing models. Again, with only a modest resolution, it will already be possible to 
make important statements. An experiment with 30° resolution could divide the emission into 
central, medial, and peripheral bins, and the relative counts would test both the concentration of 
radioisotopes in the mantle and crust, as well as the possible presence of ^'^K in the core. 

C. Uncertainties 

We want to investigate the impact of uncertainties of the reference model on the geoneutrino 
distribution. The uncertainties for crust and mantle are independent of each other. The crust itself 
can vary by about a factor of 2 in radioisotope abundances and thus in geoneutrino intensity [10]. 
The net effect of these variations thus depends on how the crust and mantle uncertainties combine. 
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FIG. 7: The plots are similar to Fig. 6, except that they show the cumulative flux, normalized as in Fig. 5. 
The dashed line is the flux of the reference model. 

Figure 8 shows some of the possible uncertainties. In the plot showing the niiniinum amount 
of radioactive elements in the crust-mantle system it can be perceived, that the reference model 
is on the low side of the possible abundances in comparison to the maximum abundances, where 
the intensity grows by a factor of two. The absence of a double peak in the intensity is due to the 
fact, that in [10] only a uncertainty for the whole crust is given, which does not take into account 
the distinction between lower, middle and upper crust. The overall shape of the maximum and 
minimum abundance plots stays overall very similar to the reference model, although there is a 
deviation from the general form at angles 9 > 30°. Thus we see that in all cases there appears a 
large peripheral flux, which remains a robust and highly testable prediction of this model. 

For a more detailed analysis it will be necessary to construct a model of the outermost layer of 
the Earth, as we assumed for a change in altitude an increase in the number of neutrinos coming 
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FIG. 8: The upper plot shows the intensity with a maximal amount of radioactive elements in crust and 
mantle. In the lower panel the intensity for a minimal abundance in crust and mantle is plotted. As the 
uncertainties in the abundances of mantle and crust are independent of each other [10], the plot in the 
middle shows a hybrid scenario with a maximal abundance in the mantle, whereas the crust abundance is 
minimized. The dashed-dotted line in all three plots is the reference model intensity is added as dashed 
dotted line in all three plots as comparison. It can be seen that the reference model is rather at the low 
side of the possible range of abundances. That the crust is only represented by a single peak is an artifact 
which arises because the uncertainties given in [10] do not distinguish between the different layers of the 
crust; this changes the intensity in the outer crust significantly. 

from the sediments. But in any case the crudeness of our estimate only infects the results for the 
outer periphery, and the central angles remain reliable, as we now see. 
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D. Crust Anisotropies and Observing Strategies 

We have assumed spherical symmetry throughout, and thus our calculations cannot directly 
address the effect of anisotropies in the radioisotope distributions. Yet these anisotropies, which 
reside in the crust, can have a very significant impact [10], since the crust is the largest radioisotope 
reservoir. For example, Mantovani et al. [10] predict fiuxes which differ by more than a factor of 2 
between locations above minimum and maximum crust depths. This would lead to changes in the 
overall geoneutrino intensity, and possibly to observable azimuthal asymmetry. 

We note that emission from the crust will affect the intensity only at the largest nadir angles. 
Thus we expect our spherical calculation to be reliable at small to median angles. Furthermore, 
one can get a rough understanding of the effect of different crust depths by adopting spherical 
models in which the crust layers have the properties appropriate for the experimental location. 

To consider an extreme example, consider a hypothetical model where we assume that the 
neutrino detector is on a ship. That means that the crust has the abundances of the oceanic crust 
and the top layer is three kilometers of ocean water. Of course in this case there is no sufficient 
shielding from cosmic rays and atmospheric radiation, which will make the detection of low energy 
neutrinos hard if not impossible. But never the less the effects are interesting, as the contribution 
coming from large angles is reduced by a factor of three, which should make it easier to detect the 
effects coming from the core and mantle. We have presently no model for an anisotropic Earth, 
but we expect, that for the real Earth in the described case the intensity for lower angles {0 < 60°) 
will be very similar to the reference model (in case of its correctness), only the contribution from 
the angles > 30° will be noticeably reduced. 

V. DISCUSSION 

Geoneutrinos offer invaluable and unique information about the energetics and structure of the 
Earth [5, 8]. The longstanding dream of measuring this neutrino population has now begun to be 
realized with the first detection of a geophysical signal by the KamLAND experiment [6]. This 
achievement is already a triumph, as the geophysical component is (by design!) dominated by the 
signal from reactor neutrinos. Nevertheless, we believe it is now worthwhile to look forward to the 
even more challenging results of determining the angular distribution of geoneutrinos, which offers 
a wealth of new information. 

In this paper we thus have calculated the angular distribution of geoneutrinos which arise 
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in /3-decays of potassium, thorium, and uranium. We have developed the general formalism for 
the neutrino intensity in a spherically symmetric Earth. We find that the geoneutrino angular 
distribution, once known, can be inverted to fully recover the terrestrial radioisotope distribution. 
Thus the geoneutrio "sky" can provide a tomography of the Earth's structure, and yields the full 
radial dependence of the radiogenic heat production. 

Turning to model-building, we explore the idealized case of an arbitrary shell of uniform density. 
This can be generalized to give an Earth model which is a series of uniform shells. We then adopt 
the radioisotope profile of the Mantovani et al. reference model [10] and calculate the resulting 
angular distribution (Fig. 4). Because the reference model places all radioisotopes in the mantle 
and crust, the resulting geoneutrino intensity is highly "peripheral," with 2/3 of the flux coming 
from nadir angles 9 > 60°. Thus, even a crude measurement of the angular distribution (say, in 
three 30° bins) would strongly test this prediction. 

We have also investigated the effect of physically plausible variations to the reference model. 
Mantovani et al. [10] identify uncertainties in their radioisotope distributions which have the effect 
of multiplicatively raising or lowering the neutrino intensity, without a significant change to the 
angular shape. However, both the shape and normalization of the intensity can change strongly 
if the Earth's core contains a significant amount of potassium, contrary to the assumptions of the 
reference model. If core potassium abundances are near the current upper limits [13] , the resulting 
geoneutrino signal can dominate the total flux. Measurement of the angular distribution will probe 
the radioisotope abundances in the core. In particular, measurement of the intensity inside a nadir 
angle 9 < 30° can discriminate among possibilities recently suggested in the literature. 

Moreover, the central and peripheral intensities are related, since the peripheral neutrinos come 
from outer shells which also contribute to the central signal. Thus a measurement of the peripheral 
flux can be used to place a lower limit to the central flux. A difference between this lower limit 
and the observed central flux then amounts to a detection of some radioisotopes in the core. 

A determination of the geoneutrino angular distribution will also solidify the connection between 
the geothermal heat flux and the geoneutrino flux. As noted by [8] , the radial component of the 
neutrino flux is directly related, by Gauss' law, to the geothermal heat production (eq. 24). But 
neutrino emission is locally isotropic and hence contains non-radial components; thus the angle- 
integrated geoneutrino flux in turn includes the non-radial components, and thus can only be related 
to the heat flux given a model of the radioisotope density profile. However, a measurement of 
the geoneutrino angular distribution can be inverted to recover the terrestrial radioisotope density 
distribution. This will allow for a full calculation not only of the global radiogenic heat production, 
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but also of its radial dependence. Thus one can test in detail models of heat production and 
transport. 

In this way, the neutrino intensity can be used to measure the radioactive contribution to the 
geothermal heat flux. This can then be compared to geophysical measurements [1] of the total 
heat flux. A comparison of these results will yield a new and robust measurement of the Urey 
ratio (eq. 1) [2]. This in turn will shed light on the thermal history of the Earth, and quantify 
the importance of any non-radioactive heating, presumably due to residual "primordial" processes 
during the formation of the Earth. 

Also, we note that the reference model we have adopted normalizes to the observed terrestrial 
heat flux and a Urey ratio of 0.5. Thus, if this Urey ratio is correct, but the heat flux contains a 
significant component from the core, this will reduce the contribution from the outer layers, which 
will act to redistribute the peripheral intensity to the interior. 

These particular results illustrate a larger more general conclusion, that geoneutrinos open a 
new window to the Earth's interior. Measurement of the angular distribution of geoneutrinos will 
allow us to infer the radioisotope distribution of the Earth. This in turn offers a new probe of 
the Earth's structure-for example, allowing a test of how sharp the radioisotope boundaries are in 
going from crust to mantle to core. If the boundaries are sharp, then the angular distribution will 
offer unprecedented new measures of the positions of these boundaries and thus will be a general 
probe of the interior structure of the Earth. 

What are the practical prospects of measuring the geoneutrino angular distribution? A complete 
discussion is beyond the scope of this paper, but it seems to us that there is some reason for 
hope. Antineutrinos are usually detected via inverse /3-decay of protons, Ug + p ^ e'^ + n, and 
the subsequent observation of the rapid e+ annihilation signal as well as the delayed n capture 
onto an ambient nucleus. Part of the momentum of the neutrino is transferred to the neutron, 
which thus contains directional information. This directionality is however lost if the neutron 
scatters before it is captured. Thus, directional sensitivity is best for liquid scintillators which use 
materials (such as gadolinium) which a large neutron capture cross section. At present, no such 
measurement appears possible: KamLAND is a liquid scintillator experiment, but the scintillator 
is an organic material which does not have a large neutron capture cross section. However, the 
reactor neutrino experiment CHOOZ did use gadolinium scintillator, and was able to reconstruct 
the source direction to within about ^ 20° [25], which as we have emphasized would already be 
geophysically interesting. While CHOOZ is no longer operating, and was too small to detect the 
geoneutrino intensity, it nevertheless demonstrated that antineutrinos with an energy spectrum 
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similar to that of geoneutrinos can be already detected with modest angular sensitivity. 

We thus strongly hope that experimental effort be undertaken to measure the angular distri- 
bution of geoneutrinos. In particular, we suggest the construction of a scintillator antineutrino 
experiment using a medium with a large neutron capture cross section, such as gadolinium. In- 
deed, there has been recent interest in somewhat similar experiments to measure the higher-energy 
signal from the cosmic supernova antineutrino background [26, 27]. While the geoneutrino flux 
is significantly higher than the supernova background, geoneutrino detection will be all the more 
challenging given that typical geoneutrino energies are relatively low, r^ 0.5 — 1.5 MeV [5, 10]. And 
of course, a good geoneutrino experiment clearly will need to be as far as possible from reactors 
(and possibly submarines!). This is particularly important if the radioisotope distribution really is 
concentrated in the outer Earth, as this will lead to a peripheral flux which can be confused with 
(but possibly calibrated by) the contribution from nearby detectors. 

We note as well that we have so far considered the total, energy-integrated, intensity. However, 
neutrino energy information is also available, and indeed geoneutrino spectra have been presented 
[10]. Since the emitted energy spectra take the well-understood /3-decay form, with sufficient energy 
resolution it is possible to separate the U, Th, and K components. If this can be done in conjunction 
with even the crudest angular resolution, it would be possible to actually distinguish beyond all 
doubt between neutrinos coming from U, Th and K and to obtain a particularly complete picture 
of the radioactive Earth. 

Thus we believe the quest to measure an image of the "geoneutrino sky" is a worthy if challenging 
goal. As we have shown, even the first, crudest attempts at neutrino imaging will yield important 
results, and will thus impel further improvements. It is our hope that this paper encourages these 
efforts. 
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APPENDIX A: TERRESTRIAL TOMOGRAPHY: INVERTING THE ANGULAR DIS- 
TRIBUTION 

We may write the intensity distribution (eq. 6) in dimensionless units as 

I{a) = 2R^ f' dx -^^= = r dx K{a, x) q{x) (Al) 

where a = sin 9 G [0, 1] and x = r/R^ G [0, 1]. Thus both / and q are defined on the interval 
[0, 1]. Furthermore, for an experiment on the surface of the Earth, we expect that /(I) = = q{l) 
because the Earth's density goes to zero at the surface (by definition!). However, a real experiment 
located slightly under the surface of the Earth might have a nonzero horizonal flux /(I). 

Clearly, I[a) is an integral transformation, with K{a, x) = xj\J x^ — cr^ the kernal. Specifically, 
eq. (Al) is a version of the Abel transform [16]. In fact, the usual Abel transform is applied to a 
function defined over an infinite domain, but fortunately one can show that the key results carry 
over to our case of a finite domain. 

The inverse Abel transform appropriate for our case is 

9(^) = -rrr / ^^ / 2 2 + „ ; 2 (^2) 

TT-K© Jx \Jo^ - x^ irRs^Vl - x^ 



^^^ d^,—IM= + ^(^ = Q) (A3) 

\/l — x"^ — iJ? ttR^\/1 — x^ 

(A4) 

where ^ = cos6, and I'{y) = dl{y)/dy is the usual derivative. 

Equation (A2) thus demonstrates by construction that, given a complete knowledge of the in- 
tensity distribution, one can fully recover the radioisotope source distribution. Thus, measurement 
of the geoneutrino angular distribution truly does carry the promise of tomographic imaging of the 
Earth's interior. In addition, with q{x) in hand, one can completely determine the radiogenic heat 
production of the Earth, both globally and as a function of depth. 

Furthermore, eq. (A2) has the properties one would expect on physical grounds. The density 
at r = Ri^x depends only on the intensity derivative for the region sin^ > x, i.e., angles along or 
exterior to the tangent angle. Thus, inferring the outer density structure requires only knowledge 
of the peripheral intensity. On the other hand, to recover the inner density structure requires both 
peripheral and central intensities. This is indeed sensible if one thinks of the angular distribution 
roughly as a linear combination of intensities along the line of sight: outer angles have only a few 
"terms" in the sum, while inner angles contain all "terms." 
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One consequence of this result is that the peripheral intensity constrains the central intensity, 
by setting a lower limit on it. If we consider I{a) only for a > a^, we can infer a lower limit Qmin 
to the density distribution at x < ctq, namely 

q{x) > qmmix) = — f da ^^ (A5) 

This example illustrates that even with an incomplete or low-resolution determination of the in- 
tensity pattern, one can draw powerful physical conclusions. 
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